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ABSTRACT 

bo ; 

We have performed a detailed V/V max test for a sample of the Canada-France 
t-h ■ Redshift Survey (CFRS) for the purpose of examining whether the comoving number 

£>v density of field galaxies changes significantly at redshifts of z < 1. Taking into 

*— h . account the luminosity evolution of galaxies which depends on their morphological 

type through different history of star formation, we obtain (V/V max ) ~ 0.5 in the 
range of 0.3 < z < 0.8, where reliable redshifts were secured by spectroscopy of either 
absorption or emission lines for the CFRS sample. This indicates that a picture of 
mild evolution of field galaxies without significant mergers is consistent with the CFRS 
(NJ ' data. Early-type galaxies, selected by their (V — I)ab color, become unnaturally 

deficient in number at z > 0.8 due to the selection bias, thereby causing a fictitious 
■ decrease of (V/V max ) . We therefore conclude that a reasonable choice of upper bound 

of redshift z ~ 0.8 in the V/V max test saves the picture of passive evolution for field 
Qs, . ellipticals in the CFRS sample, which was rejected by Kauffman, Chariot, & White 

(1996) without confining the redshift range. However, about 10% of the CFRS sample 
consists of galaxies having colors much bluer than predicted for irregular galaxies, and 
their (V/V max ) is significantly larger than 0.5. We discuss this population of extremely 

S-H ' 

blue galaxies in terms of starburst that has just turned on at their observed redshifts. 

> ' 

Subject headings: galaxies: elliptical and lenticular, cD — galaxies: evolution - 
galaxies: formation 

C3 ' 



1. Introduction 

The standard scenario of formation of elliptical galaxies is an initial starburst in dissipative 
gas collapse at very high redshift, followed by passive luminosity evolution to the present (Larson 
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1974; Tinsley & Gunn 1976). A galaxy evolution model of stellar population synthesis based on 
such a scenario well reproduces the present spectral energy distribution (SED) of elliptical galaxies 
and naturally explains their color-magnitude relation owing to the galactic wind which stops the 
starburst (Arimoto & Yoshii 1987; Matteucci & Tornambe 1987). An obvious consequence of the 
standard scenario is that the comoving number density of elliptical galaxies does not change at 
redshifts of z < 1. However, Kauffmann, Chariot, and White (1996, hereafter KCW) recently 
performed a V/V max test (Schmidt 1968) for a sample of ellipticals in the Canada-France Redshift 
Survey (CFRS, Lilly et al. 1995a), and reported a striking result that their number density 
should significantly decrease towards z ~ 1. Such a number evolution of field ellipticals obviously 
contradicts with the standard scenario and might be explained alternatively by mergers of smaller 
stellar systems and/or gaseous disks (Toomre & Toomre 1972; Kauffmann, White, & Guiderdoni 
1993; Baugh, Cole, & Frenk 1996). 

In KCWs analysis, elliptical galaxies are selected from the CFRS sample if their observed 
(V — I)ab colors are redder than a color boundary which separates ellipticals from other 
morphological types in the (V — I)ab — z diagram. They placed the boundary using a 
Bruzual-Charlot model of population synthesis with 0.1 Gyr single starburst and 50 % solar 
metallicity. This boundary certainly gives a reasonable fraction of ellipticals which agrees with 
local galaxy surveys. However, the choice of burst duration and metallicity is only ad hoc without 
any physical basis, therefore rendering some doubts as to whether the KCWs color boundary 
works successfully up to z ~ 1. 

On the other hand, local galaxies are known to have different colors from earlier to later types 
along the Hubble sequence, and this color difference is attributed to type-dependent variation 
of star formation history in galaxies (Tinsley 1980; Kennicutt, Tamblyn, & Congdon 1994). 
Therefore, the evolution models which reproduce the present colors would give a more natural way 
to distinguish various types in the (V — I)ab — z diagram. In this Letter, using the type-dependent 
evolution models of galaxies developed by Arimoto & Yoshii (1987, hereafter AY) and Arimoto, 
Yoshii, & Takahara (1992, hereafter AYT), we perform a V/V max test for the whole CFRS sample 
consisting of various types as a mixture. 

2. Color-Redshift Diagram 

Figure || shows the observed (V — I)ab colors versus secure redshifts of 591 galaxies in 
the CFRS sample. Also shown are the color trajectories based on the AY (1O 12 M0) model 
for E/SO (solid line) and the AYT (II) models for Sab (short-dashed), Sbc (long-dashed), Scd 
(dot-short-dashed), and Sdm (dot-long-dashed). Here, the symbol II for spirals stands for the 
infall model with n = 1, where n is a power index of the Schmidt law (see AYT for details). For 
comparison, KCWs model used as the color boundary between E/SO and other types is shown by 
the dotted line in this figure. 
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Unless otherwise stated, we use the A-dominated universe with (h,Q,M,^A) = (0.7,0.2,0.8), 
motivated by the fact that the luminosity density evolution based on the AY and AYT models is 
consistent with the CFRS data only in this universe (Totani, Yoshii, & Sato 1997, hereafter TYS). 
For the redshift of galaxy formation we adopt zp = 5, because the II models with this value of zp 
well reproduce Madau et al.'s (1996, 1997) observation of the UV luminosity density at z ~ 2 — 4, 
as shown by TYS and Totani (1997). Furthermore, as far as zp > 3, the effect of changing zp is 
very small in interpreting the CFRS data which extend only to z ~ 1.3. 

We see from Fig.l that the color trajectories well cover the distribution of the CFRS data, 
and this indicates an overall validity of the E/SO-Sdm models used here. In order to incorporate 
the effect of luminosity evolution in the V/V max test, we associate each of the 591 galaxies in the 
CFRS sample with a luminosity evolution obtained by interpolating between adjacent two of the 
five E/SO-Sdm models which have the nearest (V — I)ab color to a sample. We simply assign the 
E/SO or Sdm model to red or blue galaxies, respectively, which locate outside the region enclosed 
by the E/SO and Sdm color trajectories. 



3. The V/V max Test 
For a galaxy of apparent magnitude m and redshift z in the sample, we define two volumes as 

rz fjy ;-min(z„,z m ax) dV 

V = —dz and V max = / — dz , (1) 

Jmax(z ; ,z min ) UZ Jmax(z ; ,z min ) O-Z 

where the comoving volume element dV/dz is a function of z which depends on the cosmological 
parameters of h, &m, and Here, z m \ n or z max is the redshift at which the galaxy in question 
having the absolute magnitude M would be observed at the bright (m m - m ) or faint (m max ) 
magnitude limit, respectively. The CFRS is limited in apparent magnitude between m m i n = 17.5 
and m max = 22.5 in the Iab band. The redshift z x for either x='min' or x='max' should satisfy 
the following equation: 

m - mx = {K(z) + E{z)}-{K{z x )+E(z x )} + h\og^\ , (2) 

dL(Zx) 

where m and z without suffix are those of observed quantities, and dp is the luminosity distance. 
The i^-correction K and evolutionary correction E are calculated using the AY and AYT models 
for the galaxy type assigned in §2. 

Although the CFRS extends to z = 1.3, we have to limit the redshift range, z\ < z < z u , over 
which the determination of redshifts should be free from biases. Following Lilly et al. (1995b), 
we here adopt (zi,z u ) =(0.3, 1) and calculate the ratio V/V max for the CFRS galaxies in this 
redshift range. In practice, if the hypothesis of mild evolution of galaxies is correct, their values of 
V/V max are expected to distribute uniformly between and 1 with the average of V/V max equal to 
0.5, i.e., (V/Vmax) =0.5, irrespective of thier type. 
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Our V/V max statistics give (V/V milx ) =0.405 (E/SO), 0.396 (Sab), 0.507 (Sbc), 0.561 (Scd), 
and 0.642 (Sdm) with their 1 a statistical error of 0.030. A clear trend of increasing (V/V max ) from 
earlier to later types is seen, as claimed by Lilly et al (1995b) from their V/V max test without 
considering the effect of galaxy evolution. The significantly small values of (V/V max ) below 0.5 for 
E/S0 and Sab also seem to be consistent with KCW's result. We here note that all the CFRS 
galaxies have been divided into five different types of roughly equal numbers according to their 
(V — I)ab colors. Our type mixture slightly differs from local surveys, but this difference hardly 
changes the result throughout this paper. 

Since the V/V max test is known to give a biased result unless the sample is complete (Schmidt 
1968; Avni & Bahcall 1980), we should be careful in interpreting the CFRS data because of 
difficulties of redshift identification for early type galaxies at z > 1. The spectral range used for 
the CFRS covers 4250-8500 A, and there are basically no features which are useful for redshift 
identification below 3727 A at rest for emission-line galaxies and 4000 A at rest for absorption-line 
galaxies. Consequently, the redshift identification is difficult at z > 1, especially for early type 
galaxies which show the absorption lines only (Hammer et al. 1995; Crampton et al. 1995). We 
clearly see from Fig. [l| that the red, early type galaxies are unnaturally sparse at z > 0.8, while 
bluer galaxies extend beyond z ~ 1. It is therefore important to examine how the V/V max test is 
affected by changing z u in the analysis. 

Figure [2] shows (V/Vmax) as a function of z u for five different types. A steady decline of 
(V/Vmax) with increasing z u is seen beyond z u ~ 0.8 for early type galaxies due solely to their 
paucity in the range of z > 0.8. In fact, the CFRS group has also pointed out that the apparent 
absence of red quiescent objects beyond z = 0.8 is likely due to selection effects (see §7 of Hammer 
et al. 1997), which is very well consistent with our results. Therefore we have to decrease z u down 
to 0.8 in the case of performing the V/V max test free from the selection bias. In low redshifts of 
z < 0.3, on the other hand, the decrease of the CCD efficiency also makes it difficult to indentify 
red absorption galaxies (Hammer et al. 1995; Crampton et al. 1995). We have performed analysis 
varing zi, and found that changing zi leads no appreciable effect on (V ~/V max ) . Therefore we use 
the fixed value of zi = 0.3 throughout. 

Our V/V max statistics with (z u z u ) = (0.3,0.8) give (V/Vmax) =0.478 (E/S0), 0.496 (Sab), 
0.518 (Sbc), 0.513 (Scd), and 0.583 (Sdm) with their la statistical error of 0.035. These 
(V/V mauX ) values except for Sdm are consistent with 0.5 within the statistical uncertainties. 
In order to see the reason why the (V/V max ) for Sdm is significantly larger than 0.5, we 
show the CFRS galaxies classified as Sdm in Fig. ||, where V/V max is plotted against 
A(V - /)[= (V - I)% - {V - I) s /™- modcl ]. This figure indicates that much bluer galaxies 
compared to the Sdm model have larger V/V max , while galaxies which have the color nearly 
equal to the Sdm model have a relatively uniform distribution in V/V max . Apparently the large 
(^V^max) for Sdm is originated from such extremely blue galaxies (EBGs) with A(V — I) < —0.1 
as shown by the encircled points in Figs. |l] and ||. The number of EBGs is 43 in the redshift range 
of < z < 1.3, and constitutes about 10 % of the 591 galaxies in the whole CFRS sample. When 
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the EBGs are removed, {V/V max } decreases from 0.583 to 0.547 for Sdm. The (V/V max ) values still 
seem to increase from earlier to later types, and this might reflect a type transformation which is 
not included in our hypothesis of mild evolution. However, we avoid any decisive conclusion about 
this possibility of type transformation, because the {V/V max } values for all the types except for 
EBGs are consistent with 0.5 within the statistical uncertainties. Thereby we conclude that the 
hypothesis of mild evolution without significant mergers is well consistent with the CFRS data 
excluding the EBGs, in the redshift range of 0.3 < z < 0.8. 

In the above analysis we have assumed the A-dominated universe with 
(ft, Om,Oa) = (0.7,0.2,0.8). Here we repeat the calculations using the Einstein-de Sitter 
(EdS) universe with (h, Qm > ^a) = (0-5,1,0) and the low-density, open universe with 
(h, Qm, ^a) = (0.6, 0.2, 0). The results for three universe models are tabulated in Table 1, provided 
that the EBGs are removed from the sample. The average (V/V max ) increases with increasing Qm 
or decreasing Q\, and (V/V max ) of elliptical galaxies is consistent with 0.5 in all the three universe 
models. The V/V max of E/S0 in the EdS universe, 0.519, should be compared to 0.410 of KCW. 
The average (V/V max ) for the composite of all types is consistent with 0.5 in the A-dominated 
universe, and our hypothesis of mild evolution favors this universe. We furthermore repeat the 
calculations using the no-evolution model of galaxies (i.e., only X-correction is included), and the 
result is shown in the table for the purpose of comparison. It is evident that galaxy evolution 
is only a small effect in the V/V max test as far as z < 1, and our results are hardly affected by 
uncertainties of the galaxy evolution model. 



4. Discussion and Conclusions 

In addition to the selection bias in the CFRS mentioned in the previous section, our use of 
AY model for E/S0 is also responsible for the conclusion presented here which is in sharp contrast 
with KCW's conclusion. The behavior of KCW's evolution model at z < 0.6 is similar to the 
AY model, except for ~ 0.5 mag blueward shift probably due to their choice of a constant value 
of lower metallicity for all stellar populations. It should be noted that an average of absolute 
magnitude for the CFRS galaxies of E/S0 type is Mj AB = —21.9, after corrected to z = 0, 
which obviously corresponds to the AY 1O 11_12 M (baryon) models with the luminosity-weighted 
average of stellar metallicities equal to 130 — 180% of the solar (see Table 3 of AY). On the other 
hand, KCW's use of much lower metallicity of 30 — 70% of the solar is equivalent to using the AY 
1O 9_1O M0 (baryon) model with Mj AB ~ —17 which corresponds to much smaller E/S0 galaxies 
not observed in the CFRS. This suggests that KCW's evolution model is inappropriate to place a 
color boundary between E/S0 and other types for the CFRS sample. 

In fact, KCW have also noticed the selection bias against early- type galaxies at redshifts close 
to unity in the CFRS sample. Since the (V — I)ab colors are available for the CFRS galaxies with 
no redshift identifications, KCW evaluated a maximum redshift for which each unidentified galaxy 
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would still lie above the KCW's curve in Fig. 1 (the dotted line) and be classified as early type. 
They performed the V/V max test including the unidentified galaxies with these maximum redshifts 
and obtained (V/V max ) = 0.451 which is still smaller than 0.5. It should be noted that KCW's 
estimte of the maximum redshifts is heavily based on the steep rise of their model at z > 0.6 in 
Fig. 1. In contrast, our models of E/SO and Sab, which are considered to be more appropriate for 
analyzing the CFRS sample, do not show such a steep rise probably because of the longer duration 
of starbursts (0.7Gyr compared to 0.1 Gyr in KCW) and hence it is difficult or even impossible 
to define the maximum redshift for the unidentified, red galaxies. Therefore, we consider that the 
estimate of the maximum redshifts is highly uncertain, and the incompleteness should be avoided 
by performing the V/V mSuX test only in the range in which spectroscopic redshifts are secured with 
confidence. 

Let us discuss more about the nature of EBGs, which are distributed widely in the redshift 
range of < z < 1 (Fig. ||). The fact that their (V/V mauX ) value is larger than 0.5 suggests that 
the EBGs, soon after observed at z, must disappear or change their color to be classified as redder 
galaxies. A likely mechanism for such sudden fading or reddening is that the EBGs are being 
observed during intensive starburst as inferred from their extremely blue colors. In either cases 
the evolution of the EBGs may not be predicted by the usual Sdm models including ours, in which 
the star formation rate is nearly constant to make the present colors as blue as local Sdm galaxies. 
However, if the star formation is suddenly stopped on the way, the colors turn considerably redder 
than predicted by the Sdm model. It is therefore tempting to examine whether such starburst 
models for the EBGs lead to (V/V max ) ~ 0.5 for the whole CFRS galaxies. 

In this letter, we have performed a detailed V/V max test for the CFRS sample taking into 
account the AY and AYT models of galaxy evolution. Our hypothesis that all field galaxies 
evolve mildly according to type-dependent star formation histories without significant mergers is 
consistent with the CFRS data in 0.3 < z < 0.8 excluding EBGs, especially in the A-dominated 
universe. Consequently, the rejetion of passive evolution of ellipticals claimed by KCW is not 
valid. The CFRS data do not allow us to discriminate between the passive evolution and the 
number evolution of ellipticals claimed by KCW. 

If we parametrize the evolution of comoving number density of E/S0 galaxies as oc (1 + z) 7 , 
then the analysis in the limited redshift range of 0.3 < z < 0.8 gives 7e/so = 0-6 ± 1.7, —0.1 ± 1.7, 
and —0.8 ± 1.7 for the EdS, open, and A-dominated universe models, respectively. The errors are 
larger than KCW's errors, because the redshift range is limited in our analysis. For all the CFRS 
galaxies excluding the EBGs, the analysis within the same redshift range gives 7 a n = 1.8 ± 0.7, 
1.1 ± 0.7, and 0.5 ± 0.7 in the same order as above. These 7 values, though not so significant as 
claimed by KCW, would impose an important constraint on the number evolution of field galaxies 
in the framework of the hierarchical clustering of galaxies in the universe. 
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Sab 


Sbc 


Scd 


Sdm 


All Types 


(0.5, 1.0, 0.0) 


on 
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on 
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on 
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off 
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Fig. 1. — The observed redshifts versus (V — I)ab colors of the galaxies in the CFRS sample 
(Lilly et al. 1995a). The model curves are the galaxy evolution models of E/SO (solid), Sab 
(short-dashed), Sbc (long-dashed), Scd (dot-short-dashed), and Sdm (dot-long-dashed). The 
A-dominated universe with (h, Qm, ^a) = (0.7, 0.2, 0.8) and the formation redshift of zf — 5 
are assumed. The dotted line is the model used by Kauffmann, Chariot, & White (1996). 
The encircled data points are the extremely blue galaxies (EBGs) whose colors are bluer 
than the Sdm curve by more than 0.1 mag (see text for detail). 



-10- 



> 
> 

to 

> 
< 



0.6 
0.5 
0.4 

0.6 
0.5 
0.4 

0.6 
0.5 
0.4 

0.6 
0.5 
0.4 

0.6 
0.5 
0.4 





: f.T\.. 1 1 1 1 1 1 1 1 


1 


1 1 
E/SO 














1 i 1 i 1 i 1 i 




i r t 1 i 


- ' 


— h— 


1 


1 1 
Sab 


1 












i % i 1 i 1 i 1 i 1 i 












" Sbc 






— 1 — 1 — = 




■J. :.L.J 1 1 1 


— i — 


— 1 — ^ 






- — 1 — i. ■ L i. ... 1 — i — 1 — i — 




~=±=: 


Scd 












~-;J 1,1,1,1, 


, 




Sdm 

, 





0.4 0.5 0.6 0.7 0.8 0.9 1 1.1 1.2 1.3 

z„ 



Fig. 2. — The average (V/V max ) for five different types of galaxies as a function of the 
upper redshift bound z u with the fixed lower bound of z\ = 0.3. The solid and dotted lines 
represent (V/V max ) and the 1 a statistical error, respectively. The A-dominated universe 
with (h, Qm, ^a) = (0.7, 0.2, 0.8) and the formation redshift of z F = 5 are assumed. If 
(V/Vmax) is 0.5 (dashed line), the comoving number density of galaxies is conserved in the 
specified range of z\ < z < z u . 
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Fig. 3. — The V/V max versus A(V — I) plot for the Sdm galaxies selected from the the CFRS 
sample by their colors, where A(V — I) = (V — I)% - (V - I) s ^- model . The trend of larger 
V/V^x for much smaller A(V — I) can be seen. The encircled data points are the galaxies 
with A(V — I) < —0.1, which we call extremely blue galaxies (EBGs). The A-dominated 
universe with (h, Om, ^a) = (0.7, 0.2, 0.8) and the formation redshift of z F = 5 are assumed. 



